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Abstract. We present an XMM-Newton observation of the classical T Tauri star BP Tau. In the XMM-Newton RGS spectrum 
the Ovii triplet is clearly detected with a very weak forbidden line indicating high plasma densities and/or a high UV flux 
environment. At the same time concurrent UV data point to a small hot spot filling factor suggesting an accretion funnel shock 
as the site of the X-ray and UV emission. Together with the X-ray data on TW Hya these new observations suggest such funnels 
to be a general feature in classical T Tauri stars. 

Key words. X-rays: stars - stars: individual: BP Tau - stars: pre-main sequence, coronae, activity - accretion 


1. Introduction 

One of the central results of stellar X-ray astronomy is the dis¬ 
covery that all “cool stars”, i.e., stars wit h outer convective 
envel opes, are surrounded by hot coronae dSchmitt & Liefkei 
I2OO4I1 . The usual interpretation of this finding is that the com¬ 
bined action of turbulence in the outer convection zones and 
the ever-present rotation leads to dynamo action with vigor¬ 
ous production of magnetic fields and ensuing activity. X-ray 
observations support this interpretation with the “onset of ac¬ 
tivity”, i.e., the rapid increase in X-ray detection rates of main 
sequence stars at spectral type » FO and the correlation between 
rotation rate and activity. 


Strong X-ray emission is found also from many young 
stars such as zero age main sequence stars or T Tauri above 
the main sequence. Both flavors of T Tauri stars, i.e., the 
weak line T Tauri stars (wTTS) without disks as well as the 
classical T Tauri stars (cTTS) with disks are X-ray sources 
( iFeigelson & Montmerlel[r999l) . The presence of a disk around 
cTTS constitutes a fundamental difference compared to a “nor¬ 
mal” star without disk, since the energy stored in the disk 
material can be released through accretion and lead to super- 
photospheric and/or time variable emi ssions; in fact, the opti cal 
variability observed in (some) cTTS dOullbring et al.lIl996l) is 
commonly interpreted in this fashion. Yet the high-energy ac¬ 
tivity observed from cTTS and wTTS is usually interpreted as 
a “scaled-up” version of solar activity. The low-resolution X- 
ray data available prior to XMM-Newton and Chandra made a 
distinction between the X-ray properties of wTTS and cTTS 
very difficult; while some differences in emission level and 
variability properties between wTTS and cTTS were indicated 
by statistical studies dStelzer et all 2n0nHstel7er & Neiihaiisert 


l200lt iHaccomio et all I 2 OO 3 I) . a characterization of the physi¬ 
cal properties of the X-ray emitting regions was impossible. 

This situation has changed with the advent of high- 
resolution grating spectroscopy with XMM-Newton and 
Chandra. A great surprise was the X-ray spectrum of th e cTTS 
TW Hya dKastner et al.[l20()^IStelzer & Schrnitll2004l) . which 
clearly showed low forbidden to intercombination (f/i) line ra¬ 
tios in the O vii and Ne ix triplets, which were interpreted as 
due to high plasma densities in the X-ray emitti ng regions. 
The extensive spectral survey bv iNess et al ] d2004h . who stud¬ 
ied high-resolution X-ray spectra of 48 coronal sources, pro¬ 
vided no star that would even come close to the low Ovii 
f/i-ratio as observed for TW Hya. Clearly, TW Hya differs 
in its coronal properties from all other stars, and it is very 
natural to ascribe this difference to the presence of an accre¬ 
tion disk around TW Hva. S pecifically. iKastner et al 1 d2002l) 
and IStelzer & Schmitt! d2004l) interpret the X-ray emission as 
arising from an accretion shock produced by matter falling 
down onto TW Hya along a magnetic funnel. There is an obvi¬ 
ous need to examine further cTTS with high-resolution spec¬ 
troscopy in order to assess whether TW Hya’s X-ray prop¬ 
erties are typical for cTTS or not. We therefore obtained a 
high-resolution X-ray spectrum of the cTTS BP Tau with 
the Reflection Grating Spectrometer (RGS) onboard XMM- 
Newton, and the purpose of this letter is to present and discuss 
the implications of this observation. 

2. Observations 

2.1. The cTTS BP Tau: Optical and X-ray properties 

An exten sive overview of the op tical properties of BP T au is 
given by iGullbring et d] (Il996l) and lErrico et ^ (I 2 OOII) . Its 
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spectral type varies from K5-K7, it rotates fast (v sin(i) = 
15.4 km/sec), but various rotation period measurements vary¬ 
ing between 6.1 to 8.3 days have been published. The cTTS 
nature of BP Tau is clearly demonstrat ed by the excess con - 


tinuum (veiling) and Balmer emission (|Bertoute^al^ 


There is debate about BP Tau’s distance, and Wichmann et alJ 
( Il998l) argue against BP Tau being an outlier from the Taurus- 
Auriga cloud as suggested b y its HIPPARCOS para llax. From 
extensive optical monitoring Iniillbring et al.l ill 9961) conclude 
that optical variability is common in BP Tau, but in charac¬ 
ter very much different from variability encountered in typical 
flare stars. Most of the observed changes are slow and smooth 
and are interpreted as the result of inhomogeneous accretion 
from the disk onto the stellar surface. This view is strongly 
supported by the circul ar polarization in the He i A 5876 emis¬ 
sion line measured by I.Iohns-Krull et alJ lll999ll . who deduce 
a mean longitudinal magnetic field of 2460 + 120 G in the 
line forming region and argue that accretion occurs preferen¬ 
tially along large-scale magnetic loops with a small filling fac¬ 
tor. X-ray emission from BP Tau at a leve l of « 10^° erg/sec 
was first reported bv lWalter & Kuhilill98ll) using the Einstein 
Observatory. In a simult aneous optical and X-ra y observations 
of BP Tau with ROSAT iGullbring et alJ ) ll997l) found no evi¬ 
dence for correlated variations between the two bands and thus 
attributed the X-ray emission to magnetically active regions. 
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Fig. 1. OM UVW1 light curve (upper panel; bin size 300 sec), 
EPIC PN X-ray light curve (medium panel), and X-ray hard¬ 
ness ratio defined as ratio between the rates in the bands 0.2-1.0 
and 1.0-10.0 keV (lower panel; both with bin size 1000 sec). 

2.2. XMM-Newton X-Ray Data 

BP Tau was observed with XMM-Newton on August 15, 2004 
for a duration of 131 ksec (obs-id 0200370101) with the RGS 
as prime instrument. The observations were performed in full- 
frame mode employing the thick filter for both the MOS and the 
PN cameras of the European Photon Imaging Camera (EPIC). 
The optical monitor (OM) was operated using the UVW 1 fil¬ 
ter with a band pass between 2500 A and 3500 A and an ef¬ 
fective wavelength of 2910 A according to the XMM-Newton 


users’ handbook dEhle et alJ (l2004l) l. All XMM-Newton data 
was analyzed with SAS version 6.0. Background conditions 
were very quiet throughout most of the BP Tau observation and 
only lOksec of data had to be screened. In Pig. [Owe show the 
EPIC PN X-ray light curve, the X-ray hardness ratio (lower 
panel), and the UVWl light curve (upper panel). An obvious 
X-ray flare with a spectral hardness increase occurred between 
ss 20 - 40ksec into the observations; the X-ray flare may be 
accompanied by a much longer lasting UV event, but there is 
no strict correlation between X-ray and UV variability. 


BP Tou (black), TW Hya (red) 



Fig. 2. EPIC PN spectrum of BP Tau (quiescent phase; lower 
black points) in comparison to the EPIC PN spectrum of TW 
Hya (upper, red/grey points) 

In Pig. lO we show the EPIC PN spectrum of BP Tau 
(lower curve) in compa rison to TW Hya (upper curve; cf. 
IStelzer & Schmitll l2004l) . Above a; 2.5 keV the two spectra 
overlap, while at lower energies the flux of TW Hya exceeds 
that of BP Tau by almost an order of magnitude. Spectral mod¬ 
eling of the EPIC PN spectrum using simple multi-temperature 
fits requires absorption column densities of 1 - 2 x 10^' cm“^ 
depending on the chose n model and consistent wi th the relation 
between Vh and Eb-v )I.Tenkins & Savagjll974l) : hot temper¬ 
atures (ai 2keV) are required in contrast to TW Hya, whose 
EPIC spectrum is domin ated by cool plasma of » 0.3 keV (cf., 
IStelzer & Schmitlll2004l) . Thus the EPIC spectra of BP Tau and 
TW Hya suggest differences between the two stars. 

However, BP Tau’s instrumental background-subtracted high 
resolution RGS2 spectrum looks rather similar to that of TW 
Hya, (cf.. Pig. 0; probing the lower-energy plasma it shows 
mostly oxygen and neon lines with - at best - weak iron lines 
as well as signs of a continuum. Because of the low signal only 
a few spectral lines can be reliably detected in our RGS spec¬ 
tra. The clearly detected lines include the Oviii Lyo- line, the 
O VII triplet, the Nex Ly^ line, and the Neix triplet at 13.5 A; 
interestingly, no unambiguous RGS detections of iron lines, in 
particular at 15.03 A and 17.07 A were obtained (cf.. Pig. Hand 
TableO, reminiscent of the low amount of iron derived for TW 
Hya. Also, no lines of nitrogen and carbon could be detected, 
presumably because of the larger absorption column towards 
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Fig. 3. Comparison of instrumental background subtracted Fig. 4. RGS spectrum BP Tau; O vii triplet region with best fit 
RGS2 spectra of TW Hya and BP Tau with line identihcations. model. 


BP Tau. The detected lines together with the derived best ht line 
counts and their errors are listed in TableHl all line counts wer e 
derived with the CORA program dNess & Wichmanrl l2002i) . 
assuming Lorentzian line shapes. The line hts were carried 
out in such a way that the wavelength differences between the 
triplet lines and line widths were held constant; the background 
was adjusted by eye, different choices of background lead to 
line flux variations well within the erTors; in Fig.|4]we plot the 
RGS spectrum around the Ovii triplet at 22 A together with 
our best ht model. From the numbers listed in Table^we hnd 
an observed f/i-ratio of 0.37 + 0.16 for O vii and 0.40 + 0.26 
for Ne ix; the latter assumes negligible contamination by iron, 
which is a seve re problem for th e interpretation of any Neix 
triplet data (cf., iNess et al l l2003l) depending on the strength of 
the iron lines. 


Line ID 

4(A) 

Counts 

Error 

Instr. 

Flux 

Ne X Lja 

12.14 

91.5 

11.7 

RGS2 

13.2 

Ne IX He r 

13.46 

37.0 

9.9 

RGS2 

5.2 

Ne IX He i 

13.56 

31.2 

9.8 

RGS2 

4.4 

Ne IX He f 

13.71 

12.5 

7.1 

RGS2 

1.7 

O VIII Ly/; 

16.03 

25.8 

6.4 

RGSl 

4.1 

O VIII Ly/; 

16.00 

16.7 

6.2 

RGS2 

2.3 

O VIII Lyt, 

18.97 

95.5 

11.4 

RGS2 

14.2 

0 VII He r 

21.6 

47.7 

8.6 

RGSl 

8.9 

O VII He i 

21.8 

36.6 

7.8 

RGSl 

7.1 

O VII He f 

22.1 

13.7 

5.5 

RGSl 

2.7 

Fe XVII 

15.03 

< 20 


RGSl 

<2.9 


Table 1. X-ray lines detected in the RGS spectrum of BP Tau, 
line counts and photon huxes in units of 10 ® ph s * cm“^ 


3. Discussion 

IStelzer & Schmittj ( l2004l) interpret the XMM-Newton X-ray 
data on the cTTS TW Hya in terms of an accretion funnel sce¬ 
nario, where the X-ray emission is emitted in a shock (“hot 
spot”) produced by the infall of material along the magnetic 
held at essentially free-fall velocity onto the stellar surface. 


Does the same scenario also apply to the cTTS BP Tau and 
possibly to cTTS as a class ? 

It is instructive to compare the observed line huxes from 
BP Tau with those measure d from TW Hya (cf., Table 2 
in Istelzer & Schmitt! l2004l) . taking into account the differ¬ 
ent amounts of interstellar absorption. For ISM columns of 
3 X 10^°cm“^ and 1 x 10^'cm^^ respectively, we compute 
ISM transmissivities of (0.83, 0.54) at 18.97 A, (0.77, 0.42) at 
21.6 A, (0.93, 0.78) at 12.14 A, and (0.91, 0.73) at 13.46A for 
TW Hya and BP Tau respectively. These values lead to ratios 
between the Lyo, and He-like r hux ratios 1.9 and 2.1 for oxygen 
and 0.3 and 2.7 for neon respectively. While the ratios for oxy¬ 
gen are about similar, those for neon differ quite signihcantly, 
thus the emission measure distribution in BP Tau must also be 
different. We also note that the Ly^ and He-like r hux ratios ob¬ 
served in BP Tau for oxygen and neon are not consistent with 
a single temperature. 


Line 

low state 

high state “low state”/3.2 

Nex Lya 

48.6 ± 8.7 

42.7 ± 8.0 

15.2 ±2.7 

O VIII Ly„ 

59.2 ±9.1 

32.5 ± 6.6 

18.5 ±2.8 

O VII He r 

34.4 ± 7.3 

10.7 ±4.1 

10.8 ± 2.3 

O VII He 1 

24.6 ± 6.3 

9.5 ± 4.2 

7.7 ± 2.0 

O VII He f 

5.5 ±3.8 

5.6 ±3.2 

1.7 ± 1.2 


Table 2. Line counts measured in “low state” and in “high 
state” and extrapolated “high state” counts. 

We next checked to what extent individual RGS-detected 
lines are inhuenced by the hare. We dehned a “high state” by 
including all data recorded in the time interval 13.5 - 43.0 ksec 
(counted from the start of the observations as in Fig.[3i and a 
“low state” including all the rest and determined line counts as 
before. We also computed the expected number of “high state” 
counts extrapolating from the “low state” values and list all re¬ 
sults in Table 12] Clearly, within the errors the Ovii line counts 
are unaffected by the haring emission, and only 15 % of the 
recorded Oviii line counts can be attributed to the Hare; for 
Ne X this contribution rises to 30 %. We thus conclude that the 
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observed O vii emissions are not significantly affected by the 
flare, while O viii Ly^, is contaminated by < 20 %. 


In addition to the X-ray data concurrent photometric UV 
data from the XMM-Newton OM are available for BP Tau 
(cf.. Fig. The UV continuum flux of BP Tau is thought 
to be produced by an accr etion hot spot as d emonstrated by 
Ardila & Basril ( I2OOOI) and lErrico et alJ 1I2OOII) . The same ap¬ 


plies of course also to our OM UVW1 photometric data, which 
would require hlling factors of ten or more stellar surfaces if 
they were to be reconciled with BP Tau’s photospheric temper¬ 
ature. We therefore assume in the following that the recorded 
UVWl flux is exclusively produced i n a hot spot and dered- 
den this UV-flux using Ay = 0.51 dOullbring et ^ 119981) 
and Fsb-v = 0.32 x Ay. Without any further spectral in¬ 
formation we assume an emergent blackbody spectrum with 
variable temperature convert count rate into specihc en¬ 
ergy flux (erg/cm^/s/A) using the instrumental conversion fac¬ 
tor of 4.4 X 10“^® (erg/cm^/count/A), and compute the frac¬ 
tional surface area and hlling factor required to account for the 
(dereddened) observed UVWl hux. Ass uming hot spot tem¬ 
peratures between 8000 - l OOOOK (cf., ICalvet & Gullbrin'3 
Il998t lArdila & Basril Eoooli . we invariably hnd hlling factors 
between 0.6 - 4%, which are in line with values computed from 
more sophist icated models including sh ocks and an irradiated 
photosphere dCalvet & Gullbrii^ll998l) . 


Converting the observed oxygen f/i-ratio of 0.37 + 0.16 to 
density yields a nominal value of log rig - \ \ .48 with a (formal) 
error range of (11.3-11.8), neglecting radiative deexcitation of 
the forbidden line level; therefore the derived densities may be 
regarded as upper limits if the X-ray emitting region is affected 
by the observed UV radiation. We re-emphasize that both BP 
Tau and TW Hya have lower f/i-ratios and high er densities 
than a ll of the coronal X-ray sources analyzed by iNess et all 
d2004l) . If the O VII f/i-ratio was contaminated by a lower den- 
sisty “coronal component”, any accretion-related shock com¬ 
ponent would require an even smaller f/ir-ratio. 


Assuming optically thin emission, the observed energy 
flux in the Ovii and Oviii lines, /oxy, is given by /oxy = 
AshZ«e^oxy/47r£/^, where A^h, z, n^. Poxy denote shock area, 
cooling zone thickness, electron density and line cooling func¬ 
tion respectively. The post shock plasma cooling time t is 
given by T = 3nekTItigP^ot, with k, T and Ptot denoting 
Boltzmann’s constant, temperature and overall radiative loss 
function. Postshock velocity Vpost, z and t are related through 
T - z/Upost- Thus we calculate the mass accretion rate M^cc = 
'”h denoting hydrogen mass and the mean 
molecular weight. With the observed values for foxy we find 
with T = 2.5 MK Mace ~ 9 X 10“'° Mo/yr, about an order of 
magnitude smaller than inferred at UV wavelengths; since the 
observed oxygen fluxes may contain non-accretion related con¬ 
tributions, this value of M^ee should be considered as an upper 
limit. Filling factors of a few percent yield densities consistent 
with the Ovii triplet and a cooling zone thickness of » 100 
km. As to the depth of the X-ray emitting region, a model in¬ 
dependent absorption estimate using the observed photon flux 
ratio in the Lyo^ and Ly^ lines of 3.1 (cf.. Tabled requires with 
temperature of log T = 6.8 an equivalent absorption column of 


» 5 X 10^' cm“^. Given the large errors on the Ly^fLya photon 
flux ratio and the internal consistency between RGS2 values 
and optical extinction and the Ah values of Nh ~ 2xl0^*cm“^ 
derived from the broad band X-ray spectra, we conclude that 
there is no real evidence for an “additional” absorption of BP 
Tau’s X-ray flux. 

4. Conclusions 

The XMM-Newton RGS spectrum clearly demonstrates that X- 
ray emitting O vii layers in the cTTS BP Tau are either at high 
density and/or immersed in a strong ultraviolet flux. Both pos¬ 
sibilities can be well explained by an X-ray emitting accretion 
shock on BP Tau. This accretion shock can, however, produce 
only the low-temperature components in the broad band X-ray 
spectra of BP Tau (and TW Hya); to explain the high temper¬ 
ature component additional processes possibly involving mag¬ 
netic interactions between the disk and the star or magnetic 
activity anchored in the photosphere are required. Also, the de¬ 
rived mass accretion rates are smaller than those inferred from 
optical and UV data, but a detailed modeling in particular tak¬ 
ing into account non-equilibrium effects is still lacking. At any 
rate, the XMM-Newton observations of BP Tau show that TW 
Hya is not “alone”. Accretion shocks at the end of magnetic 
funnels connecting disk and stellar surface may in fact be a 
common feature of cTTS stars as a class. 
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